Broken parity and a chiral ground state in the frustrated magnet CdCr 2 04 
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We present a model describing the lattice distortion and incommensurate magnetic order in spinel 
CdCr2C>4, a good realization of the Heisenberg "pyrochlore" antiferromagnet. The magnetic frus- 
tration is relieved through the spin-Peierls distortion of the lattice involving a phonon doublet with 
odd parity. The distortion stabilizes a collinear magnetic order with the propagation wavevector 
q = 27r(0, 0, 1). The lack of inversion symmetry makes the crystal structure chiral. The handedness 
is transferred to magnetic order by the relativistic spin-orbit coupling: the collinear state is twisted 
into a long spiral with the spins in the ac plane and q shifted to 2ir(0, 5, 1). 



Frustration, denned as the presence of competing inter- 
actions, often leads to unusual effects in magnets, partic- 
ularly when it is combined with a high symmetry Q]. A 
case in point is the antiferromagnet on the "pyrochlore" 
lattice that has a very high degeneracy of the ground 
state if magnetic interactions are restricted to nearest 
neighbors |2| • At the lowest temperatures such a magnet 
is expected to retain a finite entropy per unit volume, 
as was indeed observed in a group of pyrochlore mag- 
nets known as "spin ice" t 3]. It is also well known that 
frustrated magnets with a large degeneracy are prone to 
lattice distortions that reduce the frustration by lower- 
ing the symmetry 0, 0. This effect was observed in 
antiferromagnetic spinel ZnCr204 0,0- Unfortunately, 
the distortion in this compound is rather intricate it 
involves at least four phonon modes with wavenumbers 
27r{i, |, h}. Magnetic order in the distorted lattice is 
even more complex: the magnetic unit cell is said to con- 
tain as many as 64 spins E|. As a result, the basic story 
of a flexible pyrochlore [Hi does not apply to ZnCr204 
and thus remains to be fully tested. Recent experimental 
characterization of another spinel CdCr'204 by Chung et 
al. presents us with an opportunity to do so. 

Spinels ACr204 with various nonmagnetic ions on the 
A sites are nearly ideal S — 3/2 Heisenberg antiferro- 
magnets with nearest-neighbor excha nge on the highly 
frustrated "pyrochlore" lattice |fj, [llj, [l2j. The size of 
the nonmagnetic ion determines the Cr-Cr distance and 
thereby the strength of exchange: 4.5 meV for Zn 1 
meV for Cd 01 > an d a fraction of a meV for Hg |l2| . 

OdCr^C^ undergoes a spin-Peierls- like lattice distor- 
tion 0, H at T c = 7.8 K |ll|. The lattice symmetry is 
lowered from cubic (Fdim) to tetragonal (exact space 
group unknown) with lattice constants a = b =/= c. The 
unit cell is elongated: (c — a)/c « 5 x 10~ 3 . In con- 
trast, the lattice is flattened, (c — a)/c « —1.5 x 10~ 3 , in 
ZnCr2 0*4. It is remarkable that distortions in two very 
similar compounds have opposite signs. It is also surpris- 
ing that the magnitude of the distortion is larger in the 
compound with weaker magnetic interactions. This hap- 
pens, apparently, because the quantity (c—a)/c measures 
the uniform part of the distortion only. There are indica- 



tions |fj, |j| that the nonuniform distortions in ZnC^CU 
are much larger than the uniform component. We will 
work under the assumption that the lattice distortion in 
CdCr204 lowers the point-group symmetry of the lattice 
but leaves the translational symmetry intact 10] . 

The spins in CdCr 2 04 remain disordered well below 
the Curie- Weiss temperature and order simultaneously 
with the distortion. Chung et al. \Ti\ interpret the 
magnetic order as an incommensurate spiral with the 
wavevector q = 27r(0, 5, 1) and the magnetization in the 
ac plane. They offer two ordered structures compati- 
ble with the magnetic Bragg peaks. In one the spins on 
every tetrahedron are nearly orthogonal (say, close to di- 
rections +x, +z, —x, and — z on some tetrahedra), in 
the other they are nearly collinear (say, two along +x 
and two along — x). Since 5 « 0.09 is small, we may 
treat it as an effect of a weak perturbation and begin our 
analysis at the commensurate point <5 = . 

The Landau theory of a deformable "pyrochlore" anti- 
ferromagnet yields a variety of magnetically ordered 
states, with the proposed orthogonal and collinear states 
among them. Thus, from the symmetry viewpoint, both 
candidate orders are plausible. However, a more "micro- 
scopic" treatment based on the actual physics of the spin- 
lattice coupling (and as we will see, first-principles to- 
tal energy calculations) invariably yields collinear ground 
states. Stabilization of orthogonal ground states requires 
fairly exotic interactions, such as a four-spin exchange 
strongly coupled to the lattice 0]. We therefore aban- 
don the orthogonal state and work with the collinear one. 

To simplify the calculations, we assume a clear sepa- 
ration of relevant energy scales. We treat the nearest- 
neighbor Heisenberg exchange as the strongest interac- 
tion; its minimization requires that the total spin of every 
tetrahedron be zero, which still leaves a high-dimensional 
continuum of ground states A weaker spin-lattice 
coupling selects from this continuum a collinear ground 
state. The weakest Dzyaloshinskii-Moriya (DM) interac- 
tion induces a slight misalignment of the spins and — in 
the presence of parity breaking — generates a spiral with 
a long period. In this paper we outline our findings post- 
poning a detailed account to a future publication 13]. 



2 




FIG. 1: (Color online) The collinear ground state stabilized 
by the q = 0, E u phonon (adapted from Fig. 6 of Ref. IT?}) . 
Tetrahedra of type I and II are flattened along the a and b 
directions, respectively. Solid (dashed) lines satisfied (frus- 
trated) bonds. Frustrated bonds form spirals of the same 
handedness (left-handed in this case). 

Starting point: commensurate collinear state. In a flex- 
ible pyrochlore antiferromagnet, the spin-lattice coupling 
is an efficient way to relieve spin frustration [J, ||| . The 
magnetoelastic coupling arises from the dependence of 
exchange on the ion displacements x a of spins: E me = 
(dJij/dx a ) (Si ■ Sj) x a . Integrating out the phonons gen- 
erates an effective biquadratic exchange — ^^-(Sj • Sj) 2 
favoring collinear ground states. Alternatively, this inter- 
action can be written in terms of the magnetoelastic force 
f = (/lj h) whose components are linear combinations of 
the bond variables Si ■ Sj transforming as an irreducible 
doublet E of the tetrahedral group F° r a single 

tetrahedron, the magnetoelastic energy is 
where J' is a derivative of the exchange with respect to 
ionic coordinates and k is the elastic constant of the vi- 
brational doublet E. The energy is lowest in a state with 
a tetragonal distortion, two weak and four strong bonds, 
and collinear spins [5j- 

When a distortion preserves the translational symme- 
try of the crystal, generalization to an infinte lattice is 
straightforward 10] . The existence of inequivalent tetra- 
hedra with two different orientations (I and II in Fig. 
adds inversion to the symmetry group enlarging it to 
Td ® I — Oh- The magnetoelastic energy of a primitive 
unit cell (four Cr ions) is 

E me = -K g \g\ 2 /A~K u \u\ 2 /^ (1) 



where g = f + f and u = f f are the even and 
odd doublets of bond variables whose coupling constants 
are K gu — J' 2 /k 9tU , where k g and k u are the elastic 
constants of the even and odd distortion doublets. For 
Kg > K u the lattice undergoes a uniform tetragonal dis- 
tortion with a = b > c; the space group is IAi/amd. For 
Kg < K u the distortion has both even and odd compo- 
nents: tetrahedra of types I and II are flattened along 
the a and b directions, respectively; the lattice is elon- 
gated overall, a = b < c; the space group, /4i22, lacks 
the inversion symmetry. In both cases the ground states 
are collinear (Figs. 5 and 6 in Ref. H(ih . 

The antiferromagnetic order on the pyrochlore lattice 
can be described by three staggered magnetizations Lj 
defined on tetrahedra of type I: Li = (S + Si — S 2 — 
S3)/4S' and so on. The vanishing of the total spin of a 
tetrahedron in a ground state, M 1 = X)j=o ^« = 0, makes 
the three Neel vectors Li orthogonal to each other and 
imposes a constraint on their lengths: ^ i=1 Lf = 1. In 
the state shown in Fig.^ L2 = L3 = and Li = ni e lq r 
with q = 27r(0, 0, 1); fix is an arbitrary unit vector. This 
is also consistent with the data if we take the com- 
mensurate limit 5 — * 0. 

Long-period spiral. Spiral magnetic order can arise 
when competing interactions destabilize a collinear 
ground state. That can happen when, e.g., the second- 
neighbor exchange is comparable to the nearest-neighbor 
one. However, further-neighbor exchanges are rather 
weak in spinels AC^C^ and we have checked that they 
do not destabilize collinear order (see below). 

Alternatively, spiral magnetic order may reflect a chi- 
ral nature of the underlying lattice. The handedness is 
transferred from the lattice to the spins by the relativistic 
spin-orbit coupling a(L • S). Cubic spinels are non-chiral: 
the space group FdZm includes inversion. However, par- 
ity is broken in the presence of the odd distortion E u . A 
chiral nature of the distorted lattice becomes evident if 
one examines the locations of frustrated bonds shown as 
dashed lines in Fig.^ they form spirals of the same hand- 
edness. Since the symmetry breaking is spontaneous, ex- 
periments should reveal both right and left-handed mag- 
netic spirals originating in different domains. 

In a Heisenberg magnet the spin-orbit interaction is 
manifested as the DM term Dij ■ [Si x Sj] Q. Elha- 
jal et al. [w| have determined the vectors Dy for the 
"pyrochlore" lattice up to a multiplicative constant. In 
a single tetrahedron, the DM term is 

Edm = -DS 2 (a- L 2 xL 3 +b- L 3 xLi+c- LixL 2 ). (2) 

Using the commensurate state as a starting point we 
parametrize the magnetic structure as 

L l (r)=e^(r)n l (r), (3) 

where nj(r) and 4>i(r) are the directions and magni- 
tudes of the staggered magnetizations. (Note that the 
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three unit vectors n,(r) are mutually orthogonal.) These 
parameters vary slowly in space. Proximity to the 
collinear state means that 2 and 03 are small, while 
- ((fy + 0|)/2. The staggered magnetizations © 
are defined for tetrahedra of type I. Tetrahedra of type 
II become slaves: their magnetic state is encoded in the 
staggered magnetizations of the four surrounding tetra- 
hedra of type I. The vanishing of the total magnetization 
of type-II tetrahedra yields the following constraint: 

= 03^3-6^1/4 = 0. (4) 

From it we infer that spatial derivatives of hi are of the 
same order as 02 and 03. The Neel magnetizations of a 
type-II tetrahedron are, to lowest orders, 

Li 1 = 2 n 2 - 0, Ai/4, L " = L* 1 = -8^/4. (5) 

Upon adding contributions from tetrahedra of both 
types and using Eq. (@J we obtain the DM energy 

Ebm — -DS 2 hi-(kxd x ni+bxdyiii~cxd z hi)/4:. (6) 

The terms linear in the spatial derivatives make a uniform 
state unstable against the formation of a spiral |l6| . The 
pitch of the spiral depends on the stiffness of the stag- 
gered magnetization, which is ordinarily determined by 
the strength of exchange. However, the large degeneracy 
of the pyrochlore antiferromagnet with nearest-neighbor 
exchange leads to a vanishing stiffness: indeed, apart 
from the constraint fllf . the direction of fii can vary ar- 
bitrarily in space. The stiffness is therefore determined 
by the magnetoelastic coupling and by weak exchange 
interactions beyond nearest neighbors. We discuss the 
magnetoelastic coupling first. 

A spiral magnetic state represents a deviation from the 
collinear structure and thus increases the magnetoelastic 
energy. On symmetry grounds, the increase should be 
quadratic in the gradients of fii and thus may yield a 
finite stiffness. For simplicity, we first consider only the 
odd distortions, effectively setting k g — 00 and K g = 
in Eq. JQ) , and discuss the influence of the even phonon 
later. The magnetoelastic energy can then be expressed 
in terms of the odd doublet as E mc = —K u u ■ Su/2, 
where u = AS 2 (0, 1) is the value in the commensurate 
ground state and 8u is a small deviation. As a result, we 
obtain the magnetoelastic energy density as a function of 
(j>i and the gradients of n.;. However, on account of the 
constraint Ijljl. 0q = &■* ■ d y &x/4. Likewise, minimization 
of the energy with respect to 02 yields 

02 = n 2 • d z hx/8. (7) 

The energy cost associated with the spiral is then 

E mc = K u (5 4 /4) [(cUx) 2 + (dyhrf + 2 (d 2 nx) 2 

-(na-ftni) 3 ]. (8) 




FIG. 2: (Color online) Spiral magnetic orders with the spins 
rotating about (a) the a axis, (b) the b axis, and (c) the c 
axis. 

The total energy of a spiral state is the sum of the 
DM energy © and the magnetoelastic energy JSJ. Its 
minimization yields three simple spiral states in which 
ni rotates about one of the principal axes staying in the 
plane perpendicular to it, e.g. rii = (0, cos 9 (x), sin 9 (x)). 
The energy density of all three states is the same, 

E = -DS 2 9' + K u S i 6' 2 /A, (9) 

where, in this case, 9' — d x 9. The pitch of the spiral is 

9' = 2nd = 2D/S 2 K U . (10) 

A spiral state with spins rotating about the a axis is 
shown in Fig. Ufa). On account of Eqs. Q and Q, we 
have 02 = 03 = 0, so that type-I tetrahedra still have 
collinear spins. The spins on type-II tetrahedra are copla- 
nar with the twist angle of 9'/4 = nS/2. The spiral mag- 
netic order produces Bragg scattering at q = 2n(— 8, 0, 1). 

Fig- Gib) shows a spiral state twisting about the b axis. 
This time tetrahedra of type I have coplanar spins with 
the twist angle 03 = 9'/ 4 = ttS/2, while tetrahedra of 
type II have collinear spins, as can be checked by using 
Eq. (JSJ. This spiral is related to the previous one by 
a lattice symmetry (the inversion and a n/2 rotation in 
the xy plane). It produces a magnetic Brag g p eak at 
q = 27r(0, 6, 1), as observed by Chung et al. The 
measured value 8 ~ 0.09 is consistent with a DM 
interaction that is weak relative to the magnetoelastic 
coupling: D/K U S 2 w 0.28. 

The third spiral solution is shown in Fig. Hfc). It has 
the wavevector q = 27r(0, 0,1 + 6); the spins are rotating 
around the c-axis. The twist angle is 02 = 8n/4. 

The degeneracy of the three spiral ground states is 
lifted when other perturbations are taken into account. 
A uniform E g distortion and the further-neighbor ex- 
changes J2 and J3 add energy terms 

K g (S 4 /4)[(d x ii 1 ) 2 + (dy^) 2 } + (4J 3 - 2J 2 )S 2 (c)A) 2 . 

(11) 

Depending on these coupling constants, the system will 
prefer the spiral states shown either in Fig. EJa) and (b) 
or in Fig. Efc). 
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Density-functional calculations. To test this theory, 
we performed density- functional calculations within the 
LSDA+U method using projector augmented-wave 
potentials as im ple mented in the Vienna ab initio Simu- 
lation Package Values for the on-site Coulomb 
and exchange parameters, U = 3 eV and J = 0.9 eV, 
were choosen as previously described [2(| • The results 
are not particularly sensitive to reasonable variations 
of U (±1 cV). First we performed full structural relax- 
ations in the 14-atom cubic unit cell, space group FdZm. 
Chromium ions were initialized with parallel spins in or- 
der to retain Oh symmetry throughout the structural re- 
laxation. The relaxed lattice constant, a=8.54 A, was 
found to be in excellent agreement with a=8.59 A mea- 
sured by Chung et al. fllj |. 

Estimates of exchange constants were obtained by 
comparing the total energy for several simple spin con- 
figurations. To prevent contamination by magnetoelastic 
terms lfl)l. the lattice structure was frozen in the ref- 
erence cubic state with a — 8.54 A. That procedure 
yielded J\ — 0.5 meV, J2 ~ meV, and J 3 = 0.15 
meV. The resulting Curie- Weiss temperature 0cw = 
-(l/3k B )S(S+l) J2i z i J i = -7° K compares well to the 
experimental values ranging from —70 to —90 K 

To quantify the magnetoelastic effects, we performed 
full structural relaxations for three spin configurations: 

(i) a collinear state with a pure E g distortion and mag- 
netic wavevector q = shown in Fig. 5 of Ref. llOt 

(ii) a coplanar state with a pure E g distortion and 
q = 27r(0,0, 1); and (iii) a collinear state with a mixed 
E u + Eg distortion and q = 27r(0, 0, 1) shown in Fig. ^ 
States (ii) and (iii) are the commensurate limits (5 — > 0) 
of the states displayed in Figs. 3(c) and (d) of Ref. [llj 
they have a doubled unit cell (28 atoms). The total en- 
ergy was lowest in the collinear state (iii), as posited 
above. The E g component of the distortion is tetragonal. 
Its value, (c — a)/c = +5.1 x 10 -3 , is remarkably close to 
the experimental one • Reduction of the energy asso- 
ciated with the structural relaxation depends on the spin 
state according to Eq. Q). From the data obtained in 
the three reference states we deduced the magnetoelastic 
constants K u — 0.15 meV > K g = 0.13 meV. Quantita- 
tively, the spin-lattice coupling is weaker than exchange, 
although not by much: K U S 2 / J\ = 0.67. 

The ab initio calculations back up the conclusions 
obtained analytically. The phonon doublet E u indeed 
turned out to be softer than the even distortion E g con- 
firming the selection of the collinear state of Fig. ^ The 
competition between the three candidate spiral states is 
decided by the ratio of the coupling constants in Eq. (JTTJ. 
Because J3 is quite large, the last term is prohibitively 
expensive and the spiral twists along either x or y, as 
indeed observed [ill ]. 

A large value of J3 may cast doubts on the applica- 
bility of Eq. I|ll|) , which treats that coupling as a small 
perturbation. However, it turns out that our conclusions 



in that regard remain valid for arbitrarily large values of 
J3. Luckily for us, the spiral states depicted in Fig.|5[a) 
and (b) minimize the third-neighbor exchange exactly for 
any value of the pitch 8 0|. The third spiral state in- 
creases the J3 term and is thus suppressed. 

To appreciate the unusual microscopic origin of the 
magnetic spiral, it is helpful to compare it to the standard 
scenario exemplified by the ferroelectric BiFeOs |2l|. In 
the latter case, the pitch of the magnetic spiral is propor- 
tional to an order parameter measuring the violation of 
parity, such as the electric dipolar moment. In contrast, 
in CdCr2 04 the strongest violation of parity comes not 
so much from the lattice distortion as from the magneti- 
cally ordered state itself: frustrated bonds form spirals of 
the same handedness (Fig. The corresponding order 
parameter — the odd doublet of bond variables u [lCM| — 
is a dimensionless quantity of order 1, which is why the 
presence of an order parameter is not immediately evi- 
dent in Eq. i|10|) . Without parity violation, the DM in- 
teraction alone would not generate a magnetic helix. 
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